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ABSTRACT 

Quasi-biennial oscillations (QBO) are frequently observed in the solar activity 
indices. However, no clear physical mechanism for the observed variations has 
been suggested so far. Here we study the stability of magnetic Rossby waves in the 
solar tachocline using the shallow water magnetohydrodynamic approximation. 
Our analysis shows that the combination of typical differential rotation and a 
toroidal magnetic field with a strength > 10 5 G triggers the instability of the 
m = 1 magnetic Rossby wave harmonic with a period of ~ 2 years. This harmonic 
is antisymmetric with respect to the equator and its period (and growth rate) 
depends on the differential rotation parameters and the magnetic field strength. 
The oscillations may cause a periodic magnetic flux emergence at the solar surface 
and consequently may lead to the observed QBO in the solar activity features. 
The period of QBO may change throughout the cycle, and from cycle to cycle, 
due to variations of the mean magnetic field and differential rotation in the 
tachocline. 

Subject headings: Sun: oscillations — Physical Data and Processes: magnetic 
fields — MHD — waves 



1. Introduction 

Apart from the well known 11-year cycle, solar activity shows quasi periodic variations 
on shorter time scales. Two different time scales have been frequently observed in many solar 
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activity indicators: several months and a few years. The oscillations wit h period of several 
months (mostly with 150 - 160 days) are known as Rieger-type periodicities jRieeer et al.lll984j: 
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The source(s) of these periodicities is still unclear. Several mechanisms have been sug- 
gested to dr ive the Rieger-type periodicities: interaction between I = 2 and / = 3 g-modes 
(jWolfll 119831 ), the timescale for st orage and/or escape of magnetic fields i n the solar con- 
vecti on zone (jlchimoto et al.lll985l ). "clock" modeled by an oblique rot ator (IBai fe Sturrock 
1991) and equatorially tr apped Rossby-type waves in the photosphere ( Lou 200ol ). Recently, 
Zaqarashvili et al.l ( 120101 ) (hereinafter Paper I) suggested that the Rieger-type periodicities 
can be caused by unstable m = 1, two-dimensional (8 — <fi surface in spherical coordinates) 
magnetic Rossby waves in the solar tachocline. They show that a combination of the typical 
differential rotation parameters and the magnetic field strength < 10 4 G in the tachocline 
favor the strong growth of one particular harmonic with period of 150-160 days. The pe- 
riodic modulation of the tachocline magnetic field due to the unstable harmonic triggers 
the periodic emergence of magnetic flux towards the surface, which leads to the observed 
periodicities in the solar act ivity. On the other hand, there is no clear mechanism for QBO 
reported in literature so far. iPataraya fc Zaqarashvili Jl995h supposed that the quasi 2-year 
impulse of shear waves can cause the 2-year periodicity of the differential rotation in the 
photosphere. However, this mechanism may work only near solar minima and cannot ex- 
plain the long standing modulation of solar activity. Therefore, the question of the source 
for QBO is widely open. 

In this letter, we show that the instability of magnetic Rossby waves in the tachocline 
could be the reason for QBO in solar activity. We consider a nonzero thickness of the 
tachocline and hence we use the shallow water magnetohydrodynamic (SWMHD) equations 
(IGilmanl |2000| ). We show that a stronger magnetic field (> 10 5 G) favors the growth of 
magnetic Rossby wave harmonic with period ~ 2 years. 
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2. Shallow water MHD equations and unstable magnetic Rossby wave 

harmonics 



In the following we use a spherical coordinate system (r, 9, <p) rotating with the solar 
equator, where r is the radial coordinate, 9 is the co-latitude and is the longitude. 



The solar differential rotation law in general is 
Q = f2 (l — s 2 cos 2 9 — S4 cos 1 



(i) 



where Qo is the equatorial angular velocity, and s%, S4 are constant parameters determined 
by observations. 

— * _ 

In the solar tachocline the magnetic field is predominantly toroidal, B = Se^j, and we 
take H = Ba(9) sin# , where Ba is in general a function of co-latitude. Then, the linear 
SWMHD equations (iGilmanl I200CH ) can be rewritten in the rotating frame (with Qq) as 
follows: 
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where Uq, ua, be and 6^ are the velocity and magnetic field perturbations, H is the tachocline 
thickness and h is its perturbation, p is the density, g is the reduced gravity and R is the 



distance from the solar center to the tachocline. Eqs. (E 



of shallow water magnetic field and velocity respectively ( iGilmanl 12000 ) . 



-(El) are the solenoidal conditions 



Fourier analysis with exp[im(cf) — ct)} and the transformation of variables p = cos 9 in 
Eqs. (J21)- (JSJ) lead to the equations 
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are the Alfven and surface gravity frequency respectively, h is normalized by Hq and 



#00 



b e (jl)y/l - p 2 _ Ue(lJ,)y/l - p 1 



(9) 



In the remaining we use a magnetic field 



B4, = B p, 
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which changes sign at the equator (IGilman &: Foxlll997l ). 



We expand H and h in infinite series of associated Legendre polynomials (ILonguet-Higgins 



1968|) 
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which satisfy the boundary conditions H = h = 0atfi = ±1 (i.e. at the solar poles). 
Then, we substitute Eq. (ITT]) into Eqs. ©-(JED and, using a recurrence relation of Legendre 
polynomials, we obtain algebraic equations as infinite series. The dispersion relation for 
the infinite number of harmonics can be obtained when the infinite determinant of the 
system is set to zero. In order to solve the determinant, we truncate the series at n = 
70, and the resulting polynomial in u is solved numerically. This gives the frequencies 
of the first 70 harmonics. The harmonics with real frequency are stable, while those with 



co mplex frequency are uns t able (see the ge n eral technique o 



Legendre polynomial expansion 



m lLonguet-HigginsI (Il968h : IWatsonl fll98lh : IGilman fc Foxl (119971 ): IZaqarashvili et all fl2010h 

and references therein). 

The typical values of equatorial angular velocity, radius and density in the tachocline 
are flo = 2.7 ■ 10 -6 s _1 , -Ro = 5 • 10 10 cm and p = 0.2 g • cm -3 respectively. We use a 
tachocline thickness H = 0.02_R — 10 9 cm. The ratio between angular and surface gravity 
frequencies e = Q^/Q 2 = QlRl/(gH ) is an important parameter in the shallow water theory, 
e C 1 means strongly stable stratification (main part of tachocline), while e ^> 1 considers 
weakly stable stratification (upper overshoot region). Here we consider the mean part of the 
tachocline and thus we take the limit e < 1. 



The observed differential rotation parameters near the solar surface s atisfy ~ 0.28 , 

which may tend to s 2 + s 4 ~ 0.26 near the upper part of the tachocline (jSchou et al.lll9981 ). 
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The solar radiative interior rotates uniformly, therefore the latitudinal differential rotation 
parameters drop to zero from the upper part of tachocline to its base. The radial dependence 
of latitudinal differential ro tation through th e tachocline is not clear, and S2, may also vary 
throughout the solar cycle (jliowe et al. 2000 ). Therefore, S2 + S4 may take any value between 



0.26 and 0. 

Figure [1] shows the real, mc r , and imaginary, mq, frequencies of all m = 1 unstable 
harmonics for different combinations of e (i.e. reduced gravity) and magnetic field strength. 
The differential rotation parameters are fixed to S2 = S4 = 0.11. The plot shows that each 
combination of the magnetic field strength, differential rotation parameters and reduced 
gravity leads to the occurrence of a particular harmonic whose growth rate is much stronger 
than that of other harmo nics. This is similar to what happens in the two-dimensional case 



(jZaqarashvili et al.ll2010l ). An increase of magnetic field strength leads to the reduction of 
the frequency of the most unstable harmonic. The unstable harmonics are mostly symmetric 
(defined by asterisks) with respect to the equator for a magnetic field strength < 10 4 G, while 
they become mostly antisymmetric (defined by circles) for a strength > 10 5 G. A magnetic 
field strength between 10 4 — 10 5 G yields unstable harmonics for both symmetries. This can 
be explained in terms of magnetic and differential rotation energies. Equipartition between 
the magnetic energy and the kinetic energy of differential rotation occurs at ~ 5 ■ 10 4 G for 
S2 = S4 = 0.11. When the magnetic field strength is smaller, then the differential rotation 
is the main energy source for instability and this obviously yields the symmetric harmonics 
as the differential rotation is symmetric around the equator. However, when the magnetic 
field is stronger, then the magnetic energy is the main source for the instability and the 
unstable harmonics are antisymmetric as the magnetic field is antisymmetric with respect 
to the equator. 

The importance of the equipartition value of the magnetic field strength is clearly seen on 
Figure [2J This figure displays the periods and growth rates (defined as mcj/f2 ) vs magnetic 
field strength. The growth rates are higher for weaker (< 10 4 G) and stronger (> 10 5 G) 
magnetic fields. However, the growth rates are much lower when the magnetic field strength 
is inside the interval 10 4 — 10 5 G. The weaker field (< 10 4 G) favors Rieger-type periodicities 
(150-160 days), while the stronger field (> 10 5 G) supports QBO. Increasing the magnetic 
field suppresses symmetric harmonics as it has been shown in Paper I. 

Figure [3] displays the period of the most unstable symmetric and antisymmetric har- 
monic vs the value of reduced gravity (i.e. on e) for a magnetic field strength of 8-10 4 G 
and the differential rotation parameters £2 = 54 = 0.11. It is seen that the oscillation period 
does not depend significantly on the reduced gravity. 

Figure H] displays the dependence of the periods of the most unstable symmetric and 
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antisymmetric harmonics on the differential rotation parameters for a magnetic field strength 
of 8-10 4 G and for e = 0.12 (corresponding to a reduced gravity of 1.5 • 10 2 cm s -2 ). The 
upper panel (circles) displays the antisymmetric harmonics and the lower panel (asterisks) 
displays the symmetric ones. The periods of unstable harmonics vs are plotted for different 
values of S2- The values of S2 vary from 0.13 (blue color) to 0.09 (red color). We can observe 
that the period of this harmonic depends on the differential rotation parameters significantly 
and takes values between 400-700 days. The period becomes shorter for stronger differential 
rotation. 



3. Discussion 



Our results show that the differential rotation and the magnetic field with a strength 
of > 10 5 G may lead to large-scale oscillations of tachocline with periods of ~ 2 years. The 
oscillation is due to the m = 1 unstable mode of magnetic Rossby waves. The magnetic 
Rossby waves are magnetohydr odynamic counterparts to usual hydrodynamic Rossby waves 
(jZaqarashvili et al.l 120071 . 120091 ) . The period and growth rate of the unstable harmonics 
depend on the magnetic field strength and the differential rotation parameters (Figures (U 
[2] and H]). The unstable harmonics with periods of ~ 2 years are antisymmetric with respect 
to the solar equator. 

The unstable magnetic Rossby waves in the tachocline can be the reason for QBO 
observed in almost all indices of the solar activ ity. Recent papers sugge st that QBO are 
not pers istent but may vary fro m cycle to cycle (jVecchio fc Carbonell2009l ) and throughout 
a cycle (jSykora fc Rybakl 120101 ) . Our analyses also suggest this behaviour as the period of 
unstable harmonics depends on magnetic field strength and differential rotation parameters, 
which may vary in time depending on phase and strength of a particular cycle. 

The antisymmetric behaviour of uns table harmonics with respect to the equator may 
explain the recent observational results of iBadalyan et al.l (120081 ). which show that QBO are 
well recognizable in the N-S asymmetry of solar activity indices. 

Our analysis suggests the reduction of growth rates of unstable harmonics when the 
magnetic field strength is inside the interval 10 4 — 10 5 G (see Figure [2]). It is clearly seen 
that the relatively weak magnetic field < 10 4 G leads to the occurrence of Rieger-type 
periodicities (see the same results in the paper I), while the field of > 10 5 G favors QBO. 
The upper overshoot region of the tachocline probably contains relatively weaker magnetic 
field comparing to the lower stable layers. Therefore, we may speculate that the Rieger- 
type periodicities are formed in the overshoot layer (this was also suggested in the paper 
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I), while QBO are formed in lower layers with strongly stable stratification. Therefore, the 
both periodicities may occur simultaneously. 

The magnetic field of 10 5 G is unstable due to the buoyancy instability, which makes 
difficult to keep it in the tachocline. On the other hand, the emergence of magnetic flux at 
obser ved latitud es requires sufficiently strong magnetic field (~ 10 5 ) below the convection 



zone (|Farj [200J) . Therefore, the storage of strong fields below the convection zone is still 



open question. 

Significant simplifications in our approach is the linear stability analysis, which only 
describes the initial phase of instability. Intense numerical simulations are needed in the 
future to study the developed stage of magnetic Rossby wave instabilities. 



4. Conclusions 

We have shown that the unstable magnetic Rossby waves in the solar tachocline could 
be responsible for the observed intermediate periodicities in solar activity. The periods and 
growth rates of unstable harmonics depend on the differential rotation parameters and the 
magnetic field strength. The unstable harmonics are either symmetric or antisymmetric 
with respect to the equator. The latitudinal differential rotation is mainly responsible for 
the growth of symmetric harmonics, while, the antisymmetric toroidal magnetic field favors 
the growth of antisymmetric harmonics. A magnetic field with a strength < 10 4 G leads 
to oscillations with shorter period (150-170 days), while a stronger magnetic field > 10 5 G 
favors oscillations with longer periods (1-2.5 yrs). Hence, ~ 2- year oscillations can be formed 
in the main part of the tachocline with stronger toroidal magnetic field and strongly stable 
stratification. The oscillations may trigger the periodic magnetic flux emergence at the solar 
surface and consequently QBO in solar activity. 
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Fig. 1. — Real, mc r , and imaginary, mcj, frequencies of all m = 1 unstable harmonics for 
different combinations of e (i.e. reduced gravity) and magnetic field strength. The differential 
rotation parameters are fixed to S2 = = 0.11 for all four panels. The reduced gravity varies 
from 1.5 -10 4 cm s~ 2 (e = 1.2 -10 -3 , upper left panel) to 1.1 • 10 2 cm s~ 2 (e = 0.167, lower right 
panel). Yellow, magenta, blue, green, dark blue and red colors correspond to magnetic field 
strengths of 2- 10 3 G, 2-10 4 G, 6-10 4 G, 8-10 4 G, 10 5 G,and2-10 5 G, respectively. Asterisks 
(circles) denote the symmetric (antisymmetric) harmonics with respect to the equator. 
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Fig. 2. — Period (upper panel) and growth rate mCi/Qo (lower panel) of the most unstable 
harmonics vs magnetic field strength. Asterisks (circles) define symmetric (antisymmetric) 
harmonics. Here e=0.12 and the differential rotation parameters are S2 = S4 = 0.11. 
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Fig. 3. — Periods of the most unstable symmetric (asterisks) and antisymmetric (circles) har- 
monics vs e for a magnetic field strength of 8-10 4 G and the differential rotation parameters 
s 2 = S4 = 0.11. 
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Fig. 4. — Periods of the most unstable symmetric (asterisks) and antisymmetric (circles) 
harmonics vs S4 for different values of S2- Dark blue, blue, green, magenta and red colors 
correspond to s 2 =0.13, 0.12, 0.11, 0.10 and 0.09 respectively. The magnetic field strength 
equals 8-10 4 G and e = 0.12. 



